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ABSTRACT 

This  paper  is  continued  investigations  of  propagation  of  disturbances  of  plasma 
concentration  in  the  topside  ionosphere.  A mathematical  modeling  of  the  nonlinear 
stage  of  instability  development  has  been  carried  out.  It  is  shown  that  in  the  region  of  a 
maximum  enhancement  of  perturbations,  the  relative  fluctuations  of  plasma  density  can 
make  up  several  tens  of  percent. 

In  [1,2,3]  these  authors  showed  the  possibility  of  an  enhancement  of  plasma  density 
wave  disturbances  during  their  downward  propagation  in  the  topside  ionosphere.  The 
physical  setting  of  the  problem  was  as  follows.  A harmonic  disturbance  of  plasma 
density  with  a typical  period  of  tens  to  hundreds  of  seconds  was  specified  on  a certain 
upper  boundary  (700-800  km).  Sources  for  such  disturbances  can  be  provided,  in 
particular,  by  different  types  of  oscillations  of  the  neutral  atmosphere,  such  as  AGW 
and  IGW,  and  effects  of  magnetospheric  origin.  The  propagation  of  this  disturbance 
along  geomagnetic  field  lines  was  considered.  It  was  shown  that  as  the  disturbances 
propagate  downward,  their  amplitude  can  increase  significantly,  so  that  a region  of 
strong  plasma  density  fluctuations  with  the  vertical  size  of  the  disturbances  on  the  order 
of  several  tens  of  kilometres  can  be  produced  at  about  500-600  km  altitudes.  Because 
the  spatial  growth  rates  of  enhancement  of  the  disturbances  were  quite  significant,  it  is 
of  interest  to  evaluate  the  range  of  validity  of  the  linear  approach  and  examine  nonlinear 
questions  arising  in  the  case  of  the  propagation  of  intense  disturbances. 

In  previous  work  the  basic  equation  describing  the  dynamics  of  small  disturbances  of 
electron  density  in  plasma  n,  for  the  case  of  the  ambipolar  motion  of  charges  along 
geomagnetic  field  lines,  was  derived  from  linearized  equations  of  motion  and  continuity 
for  electron-ion  gas  in  conditions  of  the  nightside  mid-latitude  ionosphere.  In  carrying 
out  the  linearization,  the  plasma  density  and  hydrodynamic  velocity  were  represented  as 
the  sum  of  the  time-independent  background  part  and  a small,  harmonically  time- 
dependent  addition:  N = N0  +n*eml ,V  = V0  +\*e'a' , and  terms  of  second  order  of 

smallness  were  n*v  discarded.  In  accordance  with  the  character  of  the  phenomenon 
under  study,  a nonlinearity  of  the  form  n*v  will  manifest  itself  in  the  generation  of 
higher  frequency  harmonics  of  external  action  co.  In  order  to  take  the  higher  harmonics 
into  account,  the  following  mathematical  description  of  the  nonlinear  processes  was 
developed.  Let  the  plasma  density  and  hydrodynamic  velocity  be  represented  as  the 
sum: 
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next,  we  substitute  (1)  into  the  equation  of  motion  and  continuity  of  electron-ion  gas: 
5V  „ 2 dN 
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dt 


= -gN  -vNW  - c 


dz 


(2) 


^L  + ^L+^N  = 0 (3) 

dt  dz 

By  combining  terms  with  identical  exponentials,  we  obtain  the  equation  for  background 
density  and  m equations  describing  harmonics  with  the  fundamental  in  m respectively. 
A natural  assumption  is  made  that  harmonic  amplitudes  decrease  with  the  harmonic 
number.  Furthermore,  higher  harmonics  are  generated  by  lower  harmonics  like  forcing 
actions.  Equations  for  the  fundamental  (4)  and  second  (5)  harmonics  are  given  below, 
and  we  may  limit  our  consideration  to  them. 
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In  the  above  equations  the  axis  z is  pointing  down,  and  the  origin  of  coordinates  lies  at 
800  km,  H is  the  height  scale  of  the  main  component  of  the  neutral  atmosphere,  atomic- 
oxygen,  Hp  is  the  plasma  height  scale,  v is  the  collision  frequency  of  ions  with  neutral 
atoms,  Vo  is  the  hydrodynamic  velocity  of  plasma  (vertical  component),  c is  the  ion 
sound  velocity,  and  p is  the  linear  recombination  coefficient.  The  quantities  v and  p 
were  assumed  to  be  exponentially  dependent  on  the  height,  and  the  velocity  Vo  was 
calculated  in  terms  of  a numerical  model  of  the  ionosphere  and  is  also  a function  of 
height.  The  other  parameters  were  taken  to  be  constant.  The  conditions  of  the  nightside 
ionosphere  of  middle  and  moderately  high  latitudes  were  considered  and  modeled. 


Figure  1. 
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Fig.  1 presents  the  fundamental  (with  the  frequency  co=0.05  c_1 ) and  second  harmonics 
in  the  linear  regime.  An  essentially  nonlinear  regime  in  a state  of  saturation  will  occur 
whenever  the  second  harmonic  becomes  comparable  with  the  fundamental  harmonic  - 
the  initial  assumption  about  decreasing  amplitudes  does  not  hold.  It  is  the  estimates  of 
this  state  that  are  of  the  greatest  interest. 

Fig.  2 plots  the  dependence  of  a maximum  of  the  ratio  of  the  second  harmonic  to  the 
fundamental  versus  ratio  of  the  fundamental  to  the  background  for  the  frequencies 
co=0.1c“'  h co— 0.01  c_1 . 


Figure  2. 

As  can  be  seen  from  the  figure,  a developed  nonlinearity  in  the  saturation  regime  will 
occur  in  the  case  of  disturbances  of  the  fundamental  harmonic  in  the  region  of  its 
maximum  on  the  order  of  20-30%  of  the  background.  Such  strong  disturbances  can 
penetrate  below  the  F2-layer  peak  and,  hence,  the  results  obtained  in  this  study  can  be 
used  to  explain  the  well-known  F-spread  phenomenon  at  vertical-incidence  radio 
soundings  of  the  mid-latitude  ionosphere. 
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